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Abstract 

Experimental data for the magnetic susceptibility and magnetostriction of 
CuGeOs are analyzed within a one-dimensional antiferromagnetic model with 
nearest (Ji) and next-nearest neighbour interactions (J2)- We show that the 
ratio of the exchange constants in the antiferromagnetic chains of CuGeO^ 
amounts to J2/J1 = 0.354(0.01), i.e. it is significantly larger than the critical 
value for the formation of a spontaneous gap in the magnetic excitation spec- 
trum without lattice dimerization. The susceptibility data are reproduced by 
our numerical results over the temperature range from 20K to 950K to a high 
degree of accuracy for Ji = 80.2(3.0) and J2 = 28.4(1.8). The pressure de- 
pendence of the exchange constants is estimated from magnetostriction data. 
Furthermore, the specific heat data are checked on consistency against the 
calculated entropy of the above model. 

PACS: 75.80.-Fq, 75.50.Ee, 75.40 Cx, 75.10.Jm 
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I. INTRODUCTION 



The investigation of low dimensional quantum spin systems has attracted widespread 
and general interest in active research on a large class of magnetic materials, both exper- 
imentally as well as theoretically, since the properties are strongly affected by quantum 
fluctuations. In particular, there is strong theoretical and experimental effort to understand 
the origin of singlet-triplet spin gaps in low dimensional spin systems such as SrCu20^ and 
V^02-P207, which occur in the presence of competing spin-interactions, e.g. due to a spin 
ladder geometry. 

Evidence for frustrated spin-interactions has also been reported for the quasi-one- 
dimensional antiferromagnet CuGeO^ HFH- Initially this compound has received much 
interest, since it is the first example of an inorganic compound undergoing a spin-Peierls 
transition |Q. Subsequent extensive experimental studies have revealed that many proper- 
ties of the ordered phase are well described by the well known spin-Peierls scenario 
such as the presence of a lattice dimerization and a singlet triplet spin-gap scaling with 
the lattice distortion However, recently it was found that the formation of the non- 
magnetic low temperature phase in CuGeO^ sensitively depends on details of the magnetic 
exchange 0. 

The traditional spin-Peierls theory is based on one-dimensional antiferromagnetic 
chains with nearest-neighbour couplings only. Such magnetic systems do not show any long- 
range antiferromagnetic order in the groundstate, however they possess critical quantum 
fluctuations which drive the lattice dimerization, i.e. the spin-Peierls transition. Substances 
with frustrated spin interactions are not strictly governed by as the magnetic system in 
these cases may show spontaneous long-range magnetic dimerization in the groundstate even 
without any lattice dimerization. The additional spin-phonon coupling merely stabilizes this 
magnetic dimerization upon developing the lattice distortion. Whether this or the former 
scenario is realized depends on the strength of the frustration parameter a = J2/ Ji- 

Frustration of the spin-interactions in CuGeOs has been inferred previously from the 
investigation of the magnetic susceptibility (x) in the non-dimerized phase , which is in 
disagreement with a nearest neighbour Heisenberg model [§,0,|T1 . A much better agreement 
has been found in theoretical studies of the spin-susceptibility x{T) invoking a Heisenberg 
chain with competing nearest and next-nearest neighbour exchange couplings Ji and J2 
ii- 

However, two markedly different choices of exchange couplings, i.e. {Ji,a) = 
(75/^,0.24) 10] and (Ji,a) = (SO/T, 0.36) were derived within model calculations and 
the same experimental data for the magnetic susceptibility as well as inelastic neutron scat- 
tering. For these quantities an increase of a on one hand and a decrease of Ji on the 
other have similar consequences leading apparently to a large uncertainty of the exchange 
parameters. 

Thus these previous studies of the quantum magnetism in CuGeO^ reveal some evidence 
for the relevance of magnetic frustration in CuGeO^, whereas it is obviously difficult to 
extract precise values of the exchange parameters. A precise knowledge of the ratio a is of 
course very important, since the theory predicts a critical ratio ac for a spin-gap to develop 
in the magnetic excitation spectrum. This gap opens irrespective of a lattice distortion. The 



existence of this gap is established exactly at the Majumdar-Ghosh point a = 1/2 [|r^] and 
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by several numerical studies [|T^,|],|l| which strongly suggest ac ~ 0.2411. 

In order to decide whether the spin gap in CuGeO^ is at least partially a manifestation of 
frustrated spin-interactions in a low dimensional magnet we have performed a comparative 
study of theoretical results and experimental data for thermodynamic properties. Using 
our numerical results for the Heisenberg chain with nearest and next-nearest neighbour 
interaction we show that a comparison of the magnetic susceptibility alone allows for the 
unique determination of the exchange constants in the frustrated one-dimensional magnet 
CuGeOs with result a = 0.354 ± 0.01, i.e. a frustration significantly larger than the critical 
value Oc- 

In the course of our investigations we also determine the pressure dependence of the 
coupling parameters from magnetostriction data. We also analyse the experimental specific 
heat data on consistency with the theoretical entropy results. 



II. THEORY AND NUMERICS 

The dominant magnetic interactions in CuGeO^ are due to Heisenberg spin exchange 
between Cu^^ ions along the c-axis of the crystal. The Hamiltonian for the spin chain in 
the non-dimerized phase reads 

H = 2 5](JiS,S,+i + J2S.S,+2), (2.1) 



where we have adopted the normalization factor of |T^. The nearest neighbour coupling Ji 



is induced by the exchange path Cu — O — Cu, and the next-nearest neighbour coupling J2 



is caused by the path Cu — O — O — Cu. A microscopic calculation of Ji 2 is difficult |[T4 
and independent derivations are important. 

In order to obtain quantitative results we calculate various physical properties in de- 
pendence on the couplings Ji and J2, notably the magnetic susceptibility, and perform a 
two-parameter fit of the experimental data which have been measured up to 950K. The aim 
is to achieve a best fit above the transition temperature Tsp = 14. 3i^ within the model 
of a magnetic system with interactions described above and an adiabatic decoupling of the 
spin-phonon interactions. In this sense Ji and J2 are treated as effective coefficients explic- 
itly dependent on the lattice geometry, i.e. microscopic bond angles and lattice constants. 
Unfortunately, analytic results for the thermodynamics of the model are available only for 
a = (nearest-neighbour Heisenberg chain). We therefore resort to complete numerical 
diagonalizations of finite systems with chain lengths up to L = 18. In general the numer- 
ical treatment of strongly correlated quantum spin chains is plagued by finite-size effects 
at low temperatures. Here, however, we are interested in relatively high temperatures with 
ksT > O.5J1. A comparison of numerical data for successive chain lengths L = 16, 17, 18 
shows that finite size corrections are essentially negligible for our purposes. In Figs. |1| 
and 1^ numerical results for the magnetic susceptibility and specific heat per lattice site 
are depicted. Note the characteristic dependence of the extremal values Xmax and the corre- 
sponding Tmax on the frustration parameter a. Tmax is decreasing with increasing a, whereas 
Xmax is increasing. The behaviour of the specific heat is similar, however its maximal value 
Cmax is a decreasing function of a. Consequently, at low temperatures the entropy strongly 
increases with increasing a. 
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We first observe that the experimental susceptibility data x allow for an unambiguous 
determination of Ji and J2 if we use the position T^ax = 56 K and absolute value x(^max) 
of the maximum of x- practice, we determine for a sequence of frustration parameters 
a = J2/J1 the value of Ji leading to Tmax = 56 K. For this set of coupling parameters the 
value of Xmax is calculated, see Fig.^]. The experimental value for Xmax is obtained in this plot 
for a = 0.354 within an error of 0.01. Note that a = 0.24 as used in would yield a value of 
Xmax far too large in comparison with the actually measured value. The nearest neighbour 
coupling corresponding to a = 0.354 is Ji = 80. 27^ ± 3.0. In Fig.^ the susceptibility x ^or 
the values a = 0.354 and 0.24 is compared with the experimental data. Note the strong 
deviation of the theoretical curve for a = 0.24 from the measured one especially at T^ax- 
In contrast to this, the entire temperature dependence of x(T) is reproduced very well for 
a = 0.354. The effect of finite-size corrections on the numerical results has been reduced 
by a scaling analysis based on a transfer matrix approach The theoretical results are 



reliable down to temperatures of 35 K. Note that a Lande factor of g = 2.256 has been 
used as derived from ESR [0. (Furthermore we have assumed a cancellation of van Vleck 



paramagnetism and core diamagnetism which is the most reasonable assumption about the 
magnetic background. If we subtract a background of 5 ■ 10~^emu/mole the result would 
read a = 0.362 ± 0.01.) For the quoted error bars of Ji, a a relative error of 3% in the 
measurements has been taken into account. Also note that a comparable analysis of the 
data m leads to a = 0.371 ± 0.01 (a = 0.38 ± 0.01 if a background of 5 ■ 10"^ emu /mole is 
subtracted), a value which overlaps with the above result. The differences between our and 
the data in ref. might originate from slightly different orientations of the magnetic fields 
with respect to the crystal axes. Note that our data agree much better with those in ref. 



In addition to the strength of the microscopic couplings we can determine their pressure 
dependence. We obtain this from magnetostriction data 0] which are related to the pressure 
dependence of the magnetic susceptibility P^PJTPj] which is rather directly accessible in 
numerical studies 

+ d,J^]], (2.2) 





TV 



where 81^2 denotes the derivative of x with respect to Ji and a. Based on this relation 
and the data for temperatures 40 K and 60 K we find dJi/dpi = 3.3(3), —7.5(5), —1.4(2) 
K GPa-\ and da/dpi = -0.03(3), 0.01(4), 0.04(2) GPa"^ for the three lattice axes i = a, 
b, and c. Note that the values for the pressure dependence of a obtained along this way 



are consistent with the hydrostatic pressure dependence obtained in [Q. However, in our 
analysis the hydrostatic pressure dependence of Ji is much stronger than that of J2 which 
is essentially zero. 



III. ENTROPY ANALYSIS 



Unfortunately it is impossible to measure the temperature dependence of the magnetic 
specific heat in CuGeOa directly. As displayed in the left part of Fig.|^ at temperatures 
close to the predicted maximum of Cmag the total specific heat is dominated by the phonon 
contribution Cph, i.e. it is about one order of magnitude larger than the calculated magnetic 
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contribution. Therefore it is impossible to extract Cmag with sufficient accuracy to resolve 
the small differences predicted for different exchange constants in a Ji — J2 model, see Fig.|^. 
Note that this would require a knowledge of the phonon contribution with an unrealistic 
accuracy of the order of 10~^ or higher. Nevertheless, one can use the measurements of C to 
check different exchange parameters suggested in the literature P JT0| , |2^ -p6| for consistency. 
In the following we will demonstrate that indeed most of these values for Ji, J2 lead to 
discrepancies. 

It is well known that a reliable separation of magnetic and phonon contributions of C 
is possible at low temperatures P^ , pB| , |?7[] and we will use this separation to estimate the 
magnetic entropy at higher temperatures. Note that the spin-Peierls phase transition does 
not modify Cph significantly, i.e. the anomaly of C is due to the magnetic contribution. This 
is inferred on one hand from the extremely small structural changes at the phase transition 
and on the other hand consistent with the findings in measurements as functions of magnetic 
fields and doping |2T]J2^. 

Analyses of the specific heat at low temperatures have been reported several times 
P^ , |23|p7| . Well below Tsp Cmag shows activated behavior due to the large spin gap. More- 
over, at low temperatures, i.e. at and below Tsp, the specific heat and correspondingly the 
entropy are dominated by their magnetic contributions. At low temperatures the phonon 
contribution follows the usual law, i.e. Cph = (3T^ with (3 ~ O.SmJ/K^mole. At higher 
temperatures Cph deviates from this behavior as illustrated in Fig.^. The extrapolation 
of the low temperature behavior exceeds the total specific heat already at moderate tem- 
peratures of about 30K. Therefore it may serve as an upper limit for Cph, which we use to 
derive a lower limit C^*^ for the magnetic contribution from the data. 

Below 20K we determine this lower limit simply from the difference C — j3T^ by assuming 
a large value of /? ~ 0.32mJ/ K'^mole. The difference function shows a maximum well above 
Tsp at 20K. A further extrapolation of this function to higher temperatures leads to a 
decrease which is unreasonable for an estimate of Cmag- We therefore take the value of 
C — PT^ at 20K as a lower limit of Cmag at temperatures above 20K (see Fig.|D. Of course 
this treatment yields only a very small lower limit for Cmag, which is indeed much smaller 
than the predicted Cmag (see Fig.||a). We do not aspire a more realistic extraction of Cmag 
here, since on one hand it is impossible to achieve the necessary accuracy to fix the exchange 
constants and on the other hand the very conservative estimate of C™^ already suffices to 
rule out certain exchange constants reported for CuGeOa. 

/^From the lower limit of Cmag we have calculated the corresponding minimum magnetic 
entropy S*,™^, which is compared to calculations of Smag for several choices of Ji and J2 in the 
right part of Fig.^. Taking the temperature at the maximum of x and assuming a = yields 
J = AAK as described in the initial paper on the spin-Peierls transition in CuGeO^ [Q. The 
corresponding entropy is much larger than the lower limit we have extracted from the data. 
Nevertheless, the specific heat data allow to exclude these parameters. At temperatures 
slightly above Tsp the calculated magnetic entropy for Ji = MK, J2 = amounts to 
about 95% of the total entropy in CuCeO^. In other words an unreasonably small phonon 
background, e.g. with a /3 ~ O.OAmJ/moleK'^ which is nearly one order of magnitude too 
small, has to be chosen to account for this large magnetic entropy. 

A further suggested value for the intrachain exchange constant in a model with a = 



is J = 60K. This value has been reported ffist in [ffUl based on their inelastic neutron 
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scattering data. More recently in ||2^ the same value and in addition a large (frustrated) 
intracliain exchange was found from the analysis of the dispersion curves in the dimerized 
phase, i.e. below Tsp. It is apparent from Fig.^ that the magnetic entropy calculated for 
these exchange parameters in a one-dimensional model, i.e. for Ji = 60K, J2 = and 
above Tsp, is significantly smaller than the lower limit S*™^ extracted from the data. Thus 



the specific heat does not support the parameters suggested in ref. |2J], in particular when 
taking into account a further reduction of the theoretical Smag due to the suggested large 
interchain exchange. 

The three solid lines in the right part of Fig.^ correspond to calculations of Smag for 
different exchange constants and finite a. Since these entropies have been calculated for 
finite chains, the results are reliable for temperatures above approx. 35K. 

The strongest discrepancy between numerical and experimental data is present for the 
set of parameters Ji = 125K, a = 0.35. These parameters have been suggested recently 
in to give the best description of the magnetic specific heat in CuGeO^,. In this latter 
work Cmag has been extracted from Raman scattering data and it was concluded that it is 
impossible to fit both the susceptibility and the magnetic specific heat with a single choice of 
Ji and J2. However, it is apparent from FigJ^ that the parameters suggested in ref. p5[ and 
consequently the "magnetic specific heat" extracted from Raman scattering is in striking 
discrepancy to the measurement of the specific heat at low temperatures. Thus our data do 
not support the reported inconsistency in the determination of exchange parameters from 



X and Cmag Further investigations seem necessary to explain the striking discrepancy 



between Cmag as revealed in pSf from quasi-elastic scattering and the true magnetic specific 
heat. 

The deviation between the data and the calculations for the exchange constants Ji = 
75i^, a = 0.24, which have been extracted by Castilla et al. from their analysis of the 
magnetic susceptibility and the dispersion curves, is less pronounced. The calculation of 
the magnetic entropy is reliable for temperatures larger than 35K where the numerical 
results are indeed larger than the lower bound we have estimated from the data. A very 
nice convergence of theory and experiment appears upon further increasing a to 0.354 (and 
J=80.2K). As shown in Fig.^ the magnetic entropy calculated for the parameters which 
yield the best fit to the susceptibility is always larger than the lower bound extracted from 
the data. Moreover, it is apparent from Fig.|^ that the difference between the theoretical 
Smag{Ji = S0.2K,a = 0.354) and the lower bound systematically decreases with decreasing 
temperature, i.e. with increasing accuracy of the extracted minimum magnetic entropy. 

Unfortunately, it is impossible to extend the theoretical calculations to temperatures of 
20K and below. At these temperatures the data for Smag shown in Fig.^ do not only represent 
a rough lower bound. Close to Tsp Smag is markedly larger than Sph for any reasonable Cph 
(As mentioned above the curve for Smag{J = 'i4:K) corresponds already to 95% of the total 
entropy). Therefore even assuming a significantly smaller phonon background leads only to 
moderate changes of Smag at these temperatures. 

Boldly extrapolating the calculated entropy for J = 80.2K, a = 0.354 to lower tem- 
peratures one may conclude that both the value of Smag at Tsp as well as its temperature 
dependence, i.e. the magnetic specific heat, are well described. In contrast to that for all 
other choices of exchange parameters in Fig.^ significant discrepancies between the model 
calculations and the experimental data are apparent. Though the accuracy for the deter- 
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mination of Cmag is not sufficient to unambiguously determine the exchange constants from 
these data alone, the specific heat strongly confirms our analysis of the susceptibility. In 
particular, in contrast to the conclusions of [^] there is no evidence that it is necessary to 
invoke markedly different exchange constants to explain x and Cmag (at low temperatures). 

IV. CONCLUSION 

We have presented numerical results for thermodynamical properties of a quantum spin- 
1/2 chain with nearest and next-nearest neighbour interactions which is believed to be 
at the heart of the magnetic system of CuGeO^. The microscopic interaction parameters 
have been determined as well as their uniaxial pressure dependence. We have shown that 
the frustration parameter is a = 0.354. This is much larger than the value used for the 
explanation of Raman scattering data We expect our result to be reliable as we 

have based our reasoning on established thermodynamical relations. Furthermore, we have 
demonstrated that the experimental magnetic susceptibility data are accounted for in even 
quantitative details by the quasi one- dimensional model and a = 0.354. 

Within the present accuracy the magnetic specific heat calculated for the exchange con- 
stants derived from our analysis of x is consistent with the analysis of the experimental 
data. On the other hand for several other choices of exchange parameters which have been 
suggested for CuGeO^ we find not only a worse description of the susceptibility but simul- 
taneously evidence for discrepancies to the specific heat data. 
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FIGURES 



FIG. 1. Plot of numerical results for the susceptibility per site (in units of Ji) versus the 
reduced temperature T/Ji for a = 0.1, 0.4, and analytical results for a = following [28| down 
to zero temperature. Shown in insets: behaviour of the susceptibility in the neighbourhood of the 
maximal values Xmax and plot of the corresponding temperature Tmax as a function of a. 



FIG. 2. Plot of numerical and analytical results for the specific heat per site similar to Fig.l. 

FIG. 3. Plot of the maximal value of the magnetic susceptibility Xmax as a function of a = J2/J1 
(solid line) in the two parameter model. The experimental value (dashed line) is crossed at 
a = 0.354. (The corresponding value of T^ax is 56K.) 

FIG. 4. Depiction of the experimental results for the magnetic susceptibility in dependence on 
temperature (solid line). Also shown are the theoretical results for a = 0.354 (dotted line) and 
a = 0.24 (dashed line). 

FIG. 5. Left upper panel: Experimentally observed specific heat of CuGeOs (o). The extrap- 
olated low temperature phonon background Cph = PT^ with f3 = 0.32mJ/ K'^mole is indicated 
by the dashed line. The solid line shows the calculated magnetic specific heat for the exchange 
constants revealing the best fit to the magnetic susceptibility. Left lower panel: Experimentally 
observed specific heat of CuGeOs (o) and the estimated minimum magnetic specific heat (solid 
line, see text). 

Right panel: Comparison between the minimum magnetic entropy as revealed from the specific 
heat data (•) and calculations of Smag assuming different exchange constants given in the figure. 
The dashed lines correspond to results of exact thermodynamic calculations for a = and the solid 
lines are obtained from numerical diagonalizations. 
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